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Abstract. Bosons in a periodic lattice with on-site disorder at low but non-zero 
temperature are considered within a mean-field theory. The criteria used for the 
definition of the superfluid, Mott insulator and Bose glass are analysed. Since the 
compressibility does never vanish at non-zero temperature, it can not be used as a 
general criterium. We show that the phases are unambiguously distinguished by the 
superfluid density and the density of states of the low-energy exitations. The phase 
diagram of the system is calculated. It is shown that even a tiny temperature leads to 
a significant shift of the boundary between the Bose glass and superfluid. 
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1. Introduction 

The remarkable experimental control over ultracold atomic gases in optical lattices 
acquired in the last couple of years [1, 2, 3, 4, 5, 6] has opened up completely new- 
lines of interest in the field of Bose-Einstein condensation. One of these are ultracold 
atoms in optical lattices with disorder which can be created by several methods. One of 
the possibilities is to use a laser speckle field [7, 8, 9]. An alternative way of creating a 
disorder potential is the introduction of a tiny fraction of a second atomic species which 
are strongly localized on random sites [10, 11, 12]. Tuning to a Feshbach resonance, these 
random scatterers can even make the disorder very strong. The random potentials 
for atoms can be also created via the spatial fluctuations of the electric currents 
generating the magnetic wire traps [13, 14, 15] or with the aid of the incommensurate 
lattices [16, 17]. Disordered lattices for ultracold rubidium atoms have been recently 
created by superimposing a regular periodic optical potential on the speckle field [18] 
and on the incommensurate lattice [19]. 

These very recent experimental developments prompt the theoretical question what 
the behavior of ultracold atomic gases in lattice potentials with disorder might be. 
Typical for the modern area of Bose-Einstein condensation (BEG) in particular and 
ultracold atomic gases in general is that such basic theoretical questions have been asked 



disordered bosons in a lattice 



2 



before in a different context. Indeed, there is an enormous body of literature on ultracold 
Fermi gases in disordered lattices referring, first of all, to electrons in amorphous solids, 
either in the normal or in the superconducting state [20]. The possibility to study 
fermionic atomic gases from this point of view is a very active and interesting field with 
which we shall not deal here, however. Suffice it to say that completely new questions 
appear like the BEC-BCS crossover [21, 22]. The influence of disorder on this problem 
is still completely unknown. 

Here, we shall focus exclusively on bosons in potentials with disorder. This is 
also known as the 'dirty boson problem'. It first came up in the pre-BEC area in 
the context of experimental investigations of superfluidity of ^He in the random pores 
of Vycor. The surprising flnding that, for sufficiently low coverage of the pores, the 
^He superfluidity would disappear even in an extrapolation to zero temperature [23, 
24, 25] prompted many theoretical studies. These were based on Hartree-Fock 
theory [26, 27], generalizations of the Bogoliubov and the Beliaev theory [28] for random 
potentials [29, 30, 31, 32, 33, 34, 35, 36], fleld-theoretical considerations [37, 38], mean- 
field theory [38, 39], renormahzation group theory [40, 41, 42, 43, 44] and quantum 
Monte Carlo simulations [45, 46, 47, 48, 49, 50, 51, 52, 53, 54] as well as numerical 
diagonalizations [55]. The consensus which developed from these studies is that in a 
disordered lattice at temperature T = two new phases of bosons may exist besides 
the supcrfluid phase which is theoretically defined by the presence of the off-diagonal 
long-range order: One is the Mott-insulator phase, which only exists at commensurate 
fillings of the lattice, and is distinguished by the absence of the off-diagonal long-range 
order, a non-zero energy gap and a vanishing compressibility. The other is the Bose- 
glass phase, which is distinguished again by the absence of the off-diagonal long-range 
order, a non-vanishing density of states at zero energy, and a non-zero compressibility. 
A more recent suggestion towards identifying the Bose-glass phase has been made in 
Ref. [56]. 

While these operational definitions of a Mott-insulator phase and a Bose-glass phase 
at temperature T = are precise and clear-cut, they run into the obvious difficulty that 
experiments and quantum Monte Carlo simulations arc never performed at T = 0. It is 
therefore necessary to examine the extent to which these or similar definitions can be 
applied at least at small non-zero temperature. This is the goal of the present paper. In 
order to achieve our goal we have to investigate the low-lying states of a suitable model 
of strongly interacting bosons in a lattice with disorder. We shall choose for this purpose 
a Bose-Hubbard model with on-site disorder of bounded variation. For our purpose of 
defining and distinguishing the various phases at non-zero temperature it is sufficient to 
analyse the basic model within a mean-field theory which allows to detect all the phases 
and to see what happens if the system parameters are changed, although it does not 
provide precise conditions for various phase transitions encountered in the system. 

The paper is organized as follows. In Section 2 the model and the mean-field 
approach to its analysis are defined. Then, in Section 3, the phase boundary between the 
superfiuid and the two non-superfluid phases is derived from the condition of vanishing 
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off-diagonal long-range order. This is first done for the pure case without disorder, 
reproducing a well-known result, and then generalizing it to the case with disorder. In 
Section 4 follows the definition of the Bose-glass and Mott-insulator phases at non-zero 
temperature and the examination of the phase boundary between them. A common 
criterion for the distinction of these two phases at T = is the compressibility. How 
well this quantity serves this purpose at finite temperature is therefore examined in 
Section 5. In Section 6 the paper ends with some final conclusions. 

2. Hamiltonian 

Wc consider a system of spinlcss bosons in a homogeneous infinitely extended lattice of 
dimension d = 1, 2, 3 described by the Bose-Hubbard Hamiltonian (in units oi h = 1) 

HbH ^-J "iOj + ^ - $Z (/^ + ^i) 4«i ' (1) 

<i,j> i i 

where J is the tunneling matrix element, U is the on-site interaction constant, and fi 
is the chemical potential. In this work we assume that the random on-site energies 
at different sites are uncorrelated and equally distributed with the probability density 

P(e). __ 

We introduce the superfluid order parameter ip — {ui), where (...) = 
Tr[. ..exp(-/3//)]/Tr[exp(-/3//)] and ... = H^ /+~ . . .p(ei)dei denote quantum- 
mechanical and disorder averaging, respectively, and f3 = l/{kT). Making use of the 
decoupling mean-field approximation in the hopping term [57, 58, 59] which is valid for 
sufficiently high-dimensional systems, we obtain the following on-site Hamiltonian 

H = -2dJ Ua^ + 'il)*a) + 2dJ + ^a^a^a a-{fjL + €) a^a , (2) 

where we have omitted the site index. 

The phase diagram of the system can be obtained in the following manner. First 
of all one has to calculate the disorder- averaged free energy of the system F{i\)) 
corresponding to the Hamiltonian (2). Then minimizing it with respect to il) to 
determine i\) = tpm one can distinguish the superfluid {ipm 7^ 0) and non-superfiuid 
(V'm = 0) regions of the parameter space. Applying a small phase gradient to the atomic 
matter-field operator and calculating the corresponding correction to the free energy in 
a manner similar to Ref. [60], one can show that the superfiuid density following from 
the Hamiltonian (2) equals [t/'ml^- In the non-superfluid region, one has to work out the 
disorder average of the static superfluid susceptibility x or the density of states p{uj) 
for the single-particle excitations [38]. In the region where p{uj) = in the interval 
Q < uj < ujg we have, by the deflnition we apply, the Mott-insulator phase with the 
energy gap Ug. On the other hand, again by deflnition, the Bose-glass phase occurs 
when lim^j^o piy^) 7^ which corresponds to a divergent superfiuid susceptibility [38]. 

The form of the mean- field Hamiltonian (2) implies that, in our approximation, 
the properties of the Mott-insulator phase as well as the Bose-glass phase, where 
vanishes, do not depend on the tunneling matrix element J. This is consistent with the 
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fact that the boundary between these phases occurs only for small values of J. The 
transition to superfluidity, where iprn starts to appear, does depend on J also in our 
approximation. 



3. Boundary between the superfiuid and non-superfiuid phases 

In order to calculate the free energy of the system, one has to solve the eigenvalue 
problem for the Hamiltonian (2). This can be done exactly by means of numerical 
calculations. However, the boundary between the superfiuid and non-superfiuid phases 
can also be determined with high accuracy treating the first term in the Hamiltonian (2) 
as a perturbation. The free energy can only depend on {ipf since a change of the phase 
in ip can be undone by the unitary transformation a ex]){—i(f)a, exp{iip)aK 
Indeed, the calculations show that the result has the following structure: 



F{^|;)^a^+a^\^|;f + ai\iJ\^ + ... (3) 

The explicit form of 04 as well as Oq and 02 for T = was obtained in Rcf. [59]. The 
generalization to T 7^ and the average over disorder needed here is straightforward. 
Since aj turns out to be always positive and can be cither positive or negative, the 
superfluid/non-supcrfiuid transition is of second order. The equation 02 = determines 
the phase boundary which is given by 

d/j,' p{n — n 



I 



H' — U{m — 1) Um — n' 



-PEmifl') ^ _J_ (A) 

2dJ ^ ' 



with 

00 jj 
Zoifx) = J2 e-^^-^'^) , EmifJ.) ^-m{m-l)-iim (5) 

m=0 ^ 

being the partition function without hopping. 



3.1. Pure case 



In the pure case we have p{e) 
2dJ g 



(5(e), so we get from (4) for the phase boundary [61] 



m=0 



m 



+ 



m + 1 
Um 



-0Em(lj) 



1 . 



(6) 



yU — f/(m — 1) um — jj, 

In the zero-tcmpcraturc limit this equation reduces to the well-known result for the 
boundary between the superfiuid and Mott-insulator [58, 59] 

\ji - U{n - 1)] [Un - ii] 



2dJ 



H + U 



(7) 



Here n denotes the positive integer at which E„i{/j) is minimal with respect to m. This 
fixes n as the smallest integer larger than or equal to fj./U. 
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3.2. Disorder with homogeneous distribution 



In the following we choose for simplicity a homogeneous disorder distribution in the 
interval e e [-A/2, A/2] 

p{e) = ^[e{e + A/2)-Q{e-A/2)], (8) 

so we have for the phase boundary (4) 



2dJ /•A'+A/2 dfj,' ~ 



m 



+ 



m+1 



ljL' — U(m — l) Um — jjt: 

We discuss first the special case T = and consider A < [/ 
ji e [U{n - 1), Un], n = 1, 2, . . .. Eq. (9) then gives 



It is assumed that 



2dJ 



nln 



II - Uin - 1) + A/2 



- (n + 1) In 



Un - It - A/2 



(10) 



// - [/(n - 1) - A/2 ■ "' Un-ii + A/2_ 

if // e [U{n - 1) + A/2, Un - A/2], otherwise 2dJ = 0. In the limit A ^ 0, Eq. (10) 
reduces to (7). The phase boundary at T = following from Eq. (10) is shown in Fig. la 
for typical parameter values. It has a lobe structure and the size of the lobes decreases 
with increasing A. If T = but A > [/, we obtain 2dJ = as the transition line for 
any value of /i, i.e., the lobes disappear and the superfluid phase appears as soon as J 
is turned on [38]. 

In the case of non-zero temperature, Eq. (9) gives always non-vanishing values 
of 2dJ. The boundary between the superfluid and non-superfluid phases for small T 
and different values of A is shown in Figs. lb,c,d. The plots obtained from Eq. (9) 
and that obtained by numerically diagonalizing the Hamiltonian (2) and minimizing 
the free energy with respect to ^0 are indistinguishable. If the temperature increases, 
the boundary between the superfluid and non-superfluid phases goes upwards and the 
size of the non-superfiuid region grows. The presence of even a small temperature 
[kT/U = 0.01) changes rather strongly the phase boundary. This is in contrast to the 
pure case where values of kT/U ~ 1 are required in order to get a noticeable shift of the 
boundary between the Mott-insulator and the superfluid phases [62]. The superfluid 
density \i/jm\ obtained by the numerical diagonalization for the same parameters as in 
Fig. lb and the values of J indicated by dotted lines in Fig. lb is plotted in Fig. 2. 



4. Boundary between the Mott-insulator and Bose-glass phases 

In the Mott-insulator as well as in the Bose-glass phase the superfluid order parameter 
ijj = iprn vanishes, which implies that J disappears from the mean-fleld Hamiltonian (2). 
Therefore, the properties of these two phases do not depend on J in the mean-field 
approximation. This is related to the fact that the two phases are localized and, 
therefore, the dependence on J should be weak. However, more accurate calculations 
beyond the mean-field theory should give some dependence on J, in particular close to 
the boundary with the superfiuid phase. The Mott-insulator phase is characterized by 
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Figure 1. Mean-field (/i, J) phase diagram for tlie homogeneous disorder 
distribution (8). At T = (a), the Bose-glass phase exists only for J = (bold lines). 
The dashed line shows the boundary (7) between the superfluid phase and the Mott- 
insulator phase in the pure case. At T 7^ (b,c,d), the Bose glass phase exists at J 7^ 
as well. If the disorder becomes strong (c,d), the Mott-insulator phase disappears. 
MI=Mott insilator with energy gap. BG=Bose glass with nonvanishing density of states 
at zero energy, SF=superfluid with nonvanishing superfluid density. The dotted lines 
in (b) indicate the two values of J used to obtain the plots in Figs. 2, 7, 8. 



0.5 
0.4 
0.3 
0.2 
0.1 



a 



0.5 



1.5 




Figure 2. Superfluid density for kT/U = 0.01, A/U = 0.5, 2dJ/U = 0.05 (a), 0.1 (b) 
(along the dotted lines in Fig. lb). 



the gap in the excitation spectrum and it has a finite superfiuid susceptibihty. The Bose- 
glass phase has no gap and the superfiuid susceptibihty diverges. AU this is directly 
related to the properties of the Green's functions and the density of states. 



4-1- Green's function 

The bosonic single-particle Green's function G{t) is defined as [63] 

G{t) = -i [e(t)G'>(t) + 0(-t)G'<(t)] , 
G>(t) = (a(t)at(0)),G<(t) = (at(0)a(t)), 



(11) 
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where a{t) is the annihilation operator in the Heisenberg representation. Straightforward 
calculations lead to the result 



-1 oo 

Zo[lJ> ) 



m=0 
oo 



G^{t) = J2 ^ e['^'-^("^-i)l^*-^^-(^') , (12) 

where fi' = n + e denotes the random local chemical potential. One can easily show 
that the imaginary-time Green's functions satisfy the periodicity condition G^{t + (3) = 
G<{t), where r = it is the imaginary time. 

At T = the imaginary-time Green's function takes the form 



(n_ + n+ + l)(n+ -n_) 1 
^>^^^ ^ 27a + VA 



{n+ + 1) e-(^-+-'^-^/2)- 



(n_ + 1) e-(^"— -hri. + n+ + 3)(n+ - n_) e"^" 



(13) 



where n± is the smallest integer greater than or equal to {ii:k A/2)/U . If n+ = n_ 
which corresponds to ji e [C/(n — 1) + A/2, C/n — A/2] for A < U, the first term in 
Eq. (13) vanishes and the superfluid susceptibihty defined by the integral over the real- 
time Green's function x — Io° G^{t)dt is a finite quantity. This means that we have not 
the Bose-glass phase, i.e., we are in the Mott-insulator phase (sec Fig. la). If n+ > n_ 
which corresponds to fi E [U{n - l),U{n - 1) + A/2] U [Un - A/2, Un] for A < f/ or 
arbitrary n for A > U, the first term in Eq. (13) survives and renders x divergent which 
is the distinguishing property of the Bose-glass phase. In the case A > U, the lobes in 
Fig. la disappear completely which means that the Mott-insulator phase is destroyed 
by the disorder [38] . 

At non-zero temperature, it is more difficult to analyze the structure of the disorder 
averaged Green's function. Expanding Eq. (12) for large but finite values of /3 shows 
that the Green's function has a similar structure as Eq. (13) but the explicit expressions 
become very long and we do not display them. Typical r-depcndcnces of in the 
Mott-insulator as well as in the Bose-glass phase arc shown for different temperatures 
in Fig. 3. Due to the different scales of r. Figs. 3a and 3b indicate that the Mott- 
insulator and the Bose-glass phase are characterized by an exponential and algebraic 
decay of G>(r), respectively. These analytical results agree qualitatively with Monte 
Carlo simulations [47, 48, 49, 52, 54]. 

Since even at finite temperature G>(t) decays only like l/r, the superfiuid 
susceptibility still diverges logarithmically. Vice versa, the divergent superfluid 
succptibility X has the consequence that the density of states for the single-particle 
excitations at zero energy does not vanish in contrast to the case of finite x [38]. In 
the next section we will see that the density of states is easier to analyze at finite 
temperature than the full Green's function. 
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(a) (b) 

Figure 3. Correlation function G> for A/U = 0.5 in the Mott- insulator [n/U = 0.5 
(a)] and Bosc-glass [ii/U = 1 (b)] phase. T = (solid hnes), kT/U = 0.01 (dashed 
lines), kT/U = 0.02 (dotted linos). Note that the scale of G'> is logarithmic but the 
scale of T is linear in (a) and logarithmic in (b), respectively. 



4-2. Density of states 



The density of states for the single-particle excitations can be determined in terms of 
the Fourier transformed single-particle Green's function G{uj) = dtexp (iut) G{t) 
as p{u,fi) = —^ImG^Lj) [28]. The Fourier transformation of Eqs. (11), (12) gives the 
density of states for the pure case 



X 



(14) 



rn5 {u + 11 — U(m — 1)) + {m + 1)5 {u + /j, — Um) 

The two 5-functions correspond to the hole and particle excitations, respectively. After 
the disorder averaging we obtain 

{m+ 1) p{Um — IX — uj) 



m=0 



ZQiUm 



^ ]^-l3Em+i{Um-uj) _|_ ^-l3Em{Um-u)) 



(15) 



This disorder averaged density of states is plotted in Fig. 4a for the Mott-insulator phase 
with the energy gap ujg = Un — A/2 — fj, and in Fig. 4b for the Bose-glass phase. Since 
Em+i{Um) — Em{Um), we get finally 



(16) 



m=o Zo{Um) 

For the homogeneous disorder distribution (8) the summation in Eq. (16) is restricted by 
m = n_, .... where n_ is the smallest integer greater than or equal to (/i — A/2)/U 
and n'_^ is the greatest integer less than or equal to (/x -|- A/2)/U. If A < t/, Eq. (16) 
takes the form 

2n ... „„„ 



-/3i5„_i(C/(n-l)) 



AZo{U{n - 1)) 





2(n+ 1) 
[ AZoiUn) 



^-l3En{Un) 



lie M 



(17) 
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Figure 4. (a) Density of states in the Mott-insulator phase for A/U = 0.5, n/U = 0.5, 
kT/U = 0.01 (solid line), 0.2 (dashed line), (b) Density of states in the Bose-glass phase 
for A/U = 0.5, fi/U = 1, kT/U = 0.01 (solid line), 0.2 (dashed line). 




Figure 5. Density of states at zero energy in the Bose-glass phase for fi/U = 1, 
A/U = 0.5 (solid line), 1 (dashed line), 1.5 (dotted line). 



where we defined 
and where Qi 



n = n{p) as the smallest integer larger than or equal p/U 
[f/(n-l),f/(n-l) + A/2], 6?2 = [f/n - A/2, f/n], M = 



[U{n — 1) + A/2, Un — A/2]. The temperature dependence of p(0, p) is plotted in Fig. 5. 
In the limit T ^ 0, Eq. (17) reduces to 

n/A p & Gi 

p{0,p) = Jo peM (18) 

^ (n+l)/A peg2. 

Simple analytical expressions (17), (18) show that the lines p = U{n — 1) + A/2 and 
p — Un — A/2 determine the boundaries between the Mott-insulator and Bose-glass 
phases at arbitrary temperature (see Fig. 1). In the case A > [/ of strong disorder, 
p(0, p) does not vanish for finite temperature and the Mott-insulator phase does not 
exist. 



5. Compressibility 

The compressibility of the system is defined as K,{p) = —d'^F{p)/dp'^, where F{p) = 
— In Z{p)/ p. In a non-superfluid phase Z{p) = Zo{p) is given by Eq. (5). Partial 
integration over e gives 

1 d r+°° dvie) 
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Figure 6. Compressibility for kT/U — 0.01, J = and the values of A/U given below 
the Unes worked out according to Eq. (20). 



For the homogeneous distribution (8) we get 



«:(/i) = -[iV(/i + A/2) 



iV (/i - A/2) 



where 



iV(/x) = (a'^a) 



— —— ^ me 



(20) 



(21) 



is the mean particle number per lattice site in the pure case. For the Bose-glass phase 
the compressibility (20) does not vanish. One can easily show that 

(3 [((ata)2) - (a^a) 

1 



lim 



J2 e-f^-^^^ - N\fi) 



If (3U ^ 1, the compressibility (20) expanded for small temperatures has the form 

-/3S 



A 



+ ae 



(22) 



(23) 



where S{fi) = En — nmi{En-i, En+i) is the energy difference between the first excited 
state and the ground state in the pure case (cf. (5)), and a is some finite constant. 
This equation shows that the Mott-insulator phase, which occurs for = n^, has an 
exponentially small compressibility at non-zero temperature, in contrast to the Bose- 
glass phase. 

The dependence of the compressibility on fi for small temperature is shown in 
Figs. 6,7. Since the compressibility does not vanish at non-zero temperature and is 
a continuous function of the system parameters, it can not be used as a criterion to 
distinguish between different phases. Thus, we deduce that the transitions are better 
defined in terms of the superfiuid order parameter ip and the density of states. 



6. Conclusions 



The Mott-insulator phase and the Bose-glass phase at vanishing temperature can be 
defined either by their thermodynamic properties or by the spectral properties of their 
low-lying excitations. Both characterizations are, of course, closely related. Both phases 
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Figure 7. Compressibility for kT/U = 0.01, A/U = 0.5, 2dJ/U = 0.05 (a), 
2dJ/U = 0.1 (b) (along the dotted lines in Fig. lb). The results are obtained by 
numerically diagonalizing the Hamiltonian (2) and calculating the second derivative of 
the free energy with respect to /x at = V'to- 



are non-supcrfluid, i.e., the corresponding Goldstone modes, the phonons, are absent. 
In the case of the Mott-insulator phase the spectral characterization by an energy gap 
imphes a vanishing compressibihty and vanishing particle-number fluctuations at T = 0. 
In the Bose-glass phase the non- vanishing density of states at zero energy implies a non- 
vanishing compressibility. These features allow a sharp distinction between the two 
phases at zero temperature. 

However, at non- vanishing temperatures, the characterization of the Mott-insulator 
and Bose-glass phases by their thermodynamic properties is no longer sharp - the Mott- 
insulator phase has an exponentially small but finite compressibility which corresponds 
to non-vanishing fiuctuations of the particle number density. Still, as we have pointed 
out in this paper, the characteristic spectral features remain present also at T > 
and can therefore be used for a sharp definition and distinction between these low 
temperature phases. We employed this possibility to calculate finite-temperature phase 
diagrams within a Bose-Hubbard model with on-site disorder within the mean-field 
approximation. For experiments with optical lattices it is usually easiest to change 
system parameters like the tunnelling amplitude J at fixed temperature, i.e., phase 
diagrams of the format of Fig. 1, where system parameters are used as variables, are 
most natural from this point of view. However, from a thermodynamic point of view, it 
may be more natural to give the phase diagram in the (yU, T)-plane. This is done in Fig. 8 
for the two values of J marked in Fig. lb as dotted lines. On the high temperature side, 
the phase diagram of Fig. 8 is incomplete, because there the transition to the normal 
gas phase must occur, which we have not considered in the present work. For both the 
Mott-insulator and the Bose-glass phase this transition would be sharp, if the energy gap 
or the finite density of states would start to appear suddenly at a critical temperature. 
Alternatively, the transition could also take the form of a smooth crossover. For the 
Mott-insulator phase the crossover could occur at the temperature kT ^ cug, where cUg 
is the energy-gap for thermal excitations [57, 62]. For the Bose-glass phase this would 
happen at the temperature where the density of states starts to be dominated by the 
normal Bose gas. 

In our calculations, we find direct transitions from the superfiuid phase either to 



disordered bosons in a lattice 

kT/U 
) . 1 1 ' 



0. 



a 



0.08 

0.06 BG 

0.04 

0.02 SF 



MI 



BG 
SF 







0.5 



12 



MI 



BG 



1.5 



0.25 
0.2 

/ 0.15 
0.1 

gp 0.05 



(b) 






MI 


BG 


BG 




BG 






SF 


Ml 

, , , 1, 


SF 





0.5 



1.5 



f,/U 



Figure 8. Phase diagram in (/x,T)-plane corresponding to the two dotted lines in 
Fig. lb, with the tunnehng ampUtude 2dJ/U = 0.05 (a), 0.1 (b). 



the Mott-insulator or to the Bose-glas phase under variation of J, depending on the 
values of // and the disorder strength A (Fig. 1). The character of the transition from 
the superfluid to the Bose-glas phase is in agreement with Monte Carlo calculations. 
The direct transition from the superfluid to the Mott-insulator phase in the case of 
weak disorder A < f/ is also in agreement with the results obtained by the path- 
integral Monte Carlo techniques [45] and by the Monte Carlo simulations based on the 
J-current model, which neglects amplitude fluctuations of the bosonic quantum fields, 
for ii/U = 0.5 [48, 49, 51, 52] and for some finite interval of /i near the tip of the Mott- 
insulator lobe [50]. However, more recent Monte Carlo investigations of the J-current 
model [53], still neglecting the amplitude fluctuations, show that the superfluid and the 
Mott-insulator regions on the phase diagram are separated by a narrow region occupied 
by the Bose-glas phase if the number of lattice sites is large enough. Thus, different 
Monte Carlo techniques give different results on the character of the superfluid - Mott- 
insulator transition in the presence of disorder. Other methods lead also to mutually 
conflicting results concerning this point [64]. Here we have presented the results of the 
mean-field approach. 

Experimentally, the Bose-glass phase may not be easy to identify with ultracold 
atoms in a suitably disordered lattice. In fact, it might best be identifiable indirectly 
by the absence of properties which are present in the competing phases for T — > 0, 
like the absence of a macroscopic wave function and the absence of an energy gap or of 
incompressibihty [19] . The finite density of states at a; ^ would show up in a specific 
heat proportional to T for T — * and in a logarithmically diverging susceptibility x- K 
would certainly be of great interest if a way could be found to measure p(0) directly. 
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